Nucleotide Exchange on RAS Proteins Using Hydrolysable and Non-hydrolysable Nucleotides.

Biochemical and biophysical assays using recombinant RAS require the protein to be in either the active or inactive state. Here we describe methods to exchange the nucleotide present in the purified RAS protein with either GDPßS, GppNHp, or GTP depending on the assay requirement. In addition, we also describe the HPLC method used to validate the exchange process and provide information on the efficiency of the nucleotide exchange.

Progress in Targeting KRAS Directly.

RAS research has entered the world of translational and clinical science. Progress has been based on our appreciation of the role of RAS mutations in different types of cancer and the effects of these mutations on the biochemical, structural, and biophysical properties of the RAS proteins themselves, particularly KRAS, on which most attention has been focused. This knowledge base, while still growing, has enabled creative chemical approaches to targeting KRAS directly. Our understanding of RAS signaling pathways in normal and cancer cells plays an important role for developing RAS inhibitors but also continues to reveal new approaches to targeting RAS through disruption of signaling complexes and downstream pathways.

Measurement of KRAS-GTPase Activity.

Oncogenic mutations in KRAS typically impact the GAP-mediated and intrinsic GTP hydrolysis activity resulting in elevated levels of cellular KRAS-GTP. The development of biochemical assays for GTPase activity provides an opportunity to quantitatively measure the impact of these mutations on GTP hydrolysis. Here we describe a biochemical assay that measures the release of free phosphate upon hydrolysis of the GTP nucleotide and allows the measurement of intrinsic or GAP-stimulated GTP hydrolysis by KRAS. This assay can be used to measure GTPase activity under single turnover conditions.

Affinity Measurement of Non-covalent Interactions of the Covalent KRAS G12C GDP Inhibitor MRTX849 to RAS Isoforms Using Surface Plasmon Resonance.

Surface plasmon resonance (SPR) is an optical effect at an electron-rich surface that enables affinity measurements of biomolecules in real time. It is label free and versatile, not limited to proteins, nucleic acids, and small molecules. SPR is a widely accepted method to measure not only affinity of molecular interactions but also association and dissociation rates of such interactions. In this chapter, we describe a general method to measure the affinity of a small molecule drug, MRTX849, to GDP bound HRAS, KRAS, and NRAS.

KRAS4b:RAF-1 Homogenous Time-Resolved Fluorescence Resonance Energy Transfer Assay for Drug Discovery.

Homogenous time-resolved FRET (HTRF) assays have become one of the most popular tools for pharmaceutical drug screening efforts over the last two decades. Large Stokes shifts and long fluorescent lifetimes of lanthanide chelates lead to robust signal to noise, as well as decreased false positive rates compared to traditional assay techniques. In this chapter, we describe an HTRF protein-protein interaction (PPI) assay for the KRAS4b G-domain in the GppNHp-bound state and the RAF-1-RBD currently used for drug screens. Application of this assay contributes to the identification of lead compounds targeting the GTP-bound active state of K-RAS.

A Top-Down Proteomic Assay to Evaluate KRAS4B-Compound Engagement.

Development of new targeted inhibitors for oncogenic KRAS mutants may benefit from insight into how a given mutation influences the accessibility of protein residues and how compounds interact with mutant or wild-type KRAS proteins. Targeted proteomic analysis, a key validation step in the KRAS inhibitor development process, typically involves both intact mass- and peptide-based methods to confirm compound localization or quantify binding. However, these methods may not always provide a clear picture of the compound binding affinity for KRAS, how specific the compound is to the target KRAS residue, and how experimental conditions may impact these factors. To address this, we have developed a novel top-down proteomic assay to evaluate in vitro KRAS4B-compound engagement while assessing relative quantitation in parallel. We present two applications to demonstrate the capabilities of our assay: maleimide-biotin labeling of a KRAS4BG12D cysteine mutant panel and treatment of three KRAS4B proteins (WT, G12C, and G13C) with small molecule compounds. Our results show the time- or concentration-dependence of KRAS4B-compound engagement in context of the intact protein molecule while directly mapping the compound binding site.

Membrane lipids drive formation of KRAS4b-RAF1 RBDCRD nanoclusters on the membrane.

The oncogene RAS, extensively studied for decades, presents persistent gaps in understanding, hindering the development of effective therapeutic strategies due to a lack of precise details on how RAS initiates MAPK signaling with RAF effector proteins at the plasma membrane. Recent advances in X-ray crystallography, cryo-EM, and super-resolution fluorescence microscopy offer structural and spatial insights, yet the molecular mechanisms involving protein-protein and protein-lipid interactions in RAS-mediated signaling require further characterization. This study utilizes single-molecule experimental techniques, nuclear magnetic resonance spectroscopy, and the computational Machine-Learned Modeling Infrastructure (MuMMI) to examine KRAS4b and RAF1 on a biologically relevant lipid bilayer. MuMMI captures long-timescale events while preserving detailed atomic descriptions, providing testable models for experimental validation. Both in vitro and computational studies reveal that RBDCRD binding alters KRAS lateral diffusion on the lipid bilayer, increasing cluster size and decreasing diffusion. RAS and membrane binding cause hydrophobic residues in the CRD region to penetrate the bilayer, stabilizing complexes through ß-strand elongation. These cooperative interactions among lipids, KRAS4b, and RAF1 are proposed as essential for forming nanoclusters, potentially a critical step in MAP kinase signal activation.

Comparative analysis of KRAS4a and KRAS4b splice variants reveals distinctive structural and functional properties.

KRAS, the most frequently mutated oncogene in human cancer, produces two isoforms, KRAS4a and KRAS4b, through alternative splicing. These isoforms differ in exon 4, which encodes the final 15 residues of the G-domain and hypervariable regions (HVRs), vital for trafficking and membrane localization. While KRAS4b has been extensively studied, KRAS4a has been largely overlooked. Our multidisciplinary study compared the structural and functional characteristics of KRAS4a and KRAS4b, revealing distinct structural properties and thermal stability. Position 151 influences KRAS4a's thermal stability, while position 153 affects binding to RAF1 CRD protein. Nuclear magnetic resonance analysis identified localized structural differences near sequence variations and provided a solution-state conformational ensemble. Notably, KRAS4a exhibits substantial transcript abundance in bile ducts, liver, and stomach, with transcript levels approaching KRAS4b in the colon and rectum. Functional disparities were observed in full-length KRAS variants, highlighting the impact of HVR variations on interaction with trafficking proteins and downstream effectors like RAF and PI3K within cells.

Revealing the mechanism of action of a first-in-class covalent inhibitor of KRASG12C (ON) and other functional properties of oncogenic KRAS by 31P NMR.

Individual oncogenic KRAS mutants confer distinct differences in biochemical properties and signaling for reasons that are not well understood. KRAS activity is closely coupled to protein dynamics and is regulated through two interconverting conformations: state 1 (inactive, effector binding deficient) and state 2 (active, effector binding enabled). Here, we use 31P NMR to delineate the differences in state 1 and state 2 populations present in WT and common KRAS oncogenic mutants (G12C, G12D, G12V, G13D, and Q61L) bound to its natural substrate GTP or a commonly used nonhydrolyzable analog GppNHp (guanosine-5'-[(ß,γ)-imido] triphosphate). Our results show that GppNHp-bound proteins exhibit significant state 1 population, whereas GTP-bound KRAS is primarily (90% or more) in state 2 conformation. This observation suggests that the predominance of state 1 shown here and in other studies is related to GppNHp and is most likely nonexistent in cells. We characterize the impact of this differential conformational equilibrium of oncogenic KRAS on RAF1 kinase effector RAS-binding domain and intrinsic hydrolysis. Through a KRAS G12C drug discovery, we have identified a novel small-molecule inhibitor, BBO-8956, which is effective against both GDP- and GTP-bound KRAS G12C. We show that binding of this inhibitor significantly perturbs state 1-state 2 equilibrium and induces an inactive state 1 conformation in GTP-bound KRAS G12C. In the presence of BBO-8956, RAF1-RAS-binding domain is unable to induce a signaling competent state 2 conformation within the ternary complex, demonstrating the mechanism of action for this novel and active-conformation inhibitor.

Revealing KRas4b topology on the membrane surface.

KRas4b is a membrane-bound regulatory protein belonging to the family of small GTPases that function as a molecular switch, facilitating signal transduction from activated membrane receptors to intracellular pathways controlling cell growth and proliferation. Oncogenic mutations locking KRas4b in the active GTP state are responsible for nearly 85% of all Ras-driven cancers. Understanding the membrane-bound state of KRas4b is crucial for designing new therapeutic approaches targeting oncogenic KRas-driven signaling pathways. Extensive research demonstrates the significant involvement of the membrane bilayer in Ras-effector interactions, with anionic lipids playing a critical role in determining protein conformations The preferred topology of KRas4b for interacting with signaling partners has been a long-time question. Computational studies suggest a membrane-proximal conformation, while other biophysical methods like neutron reflectivity propose a membrane-distal conformation. To address these gaps, we employed FRET measurements to investigate the conformation of KRas4b. Using fully post-translationally modified KRas4b, we designed a Nanodisc based FRET assay to study KRas4b-membrane interactions. We suggest an extended conformation of KRas4b relative to the membrane surface. Measurement of FRET donor - acceptor distances reveal that a negatively charged membrane surface weakly favors closer association with the membrane surface. Our findings provide insights into the role of anionic lipids in determining the dynamic conformations of KRas4b and shed light on the predominant conformation of its topology on lipid headgroups.

Direct Modulators of K-Ras-Membrane Interactions.

Protein-membrane interactions (PMIs) are ubiquitous in cellular signaling. Initial steps of signal transduction cascades often rely on transient and dynamic interactions with the inner plasma membrane leaflet to populate and regulate signaling hotspots. Methods to target and modulate these interactions could yield attractive tool compounds and drug candidates. Here, we demonstrate that the conjugation of a medium-chain lipid tail to the covalent K-Ras(G12C) binder MRTX849 at a solvent-exposed site enables such direct modulation of PMIs. The conjugated lipid tail interacts with the tethered membrane and changes the relative membrane orientation and conformation of K-Ras(G12C), as shown by molecular dynamics (MD) simulation-supported NMR studies. In cells, this PMI modulation restricts the lateral mobility of K-Ras(G12C) and disrupts nanoclusters. The described strategy could be broadly applicable to selectively modulate transient PMIs.

Allosteric Regulation of Switch-II Domain Controls KRAS Oncogenicity.

RAS proteins are GTPases that regulate a wide range of cellular processes. RAS activity is dependent on its nucleotide-binding status, which is modulated by guanine nucleotide exchange factors (GEF) and GTPase-activating proteins (GAP). KRAS can be acetylated at lysine 104 (K104), and an acetylation-mimetic mutation of K104 to glutamine (K104Q) attenuates the in vitro-transforming capacity of oncogenic KRAS by interrupting GEF-induced nucleotide exchange. To assess the effect of this mutation in vivo, we used CRISPR-Cas9 to generate mouse models carrying the K104Q point mutation in wild-type and conditional KrasLSL-G12D alleles. Homozygous animals for K104Q were viable, fertile, and arose at the expected Mendelian frequency, indicating that K104Q is not a complete loss-of-function mutation. Consistent with our previous findings from in vitro studies, however, the oncogenic activity of KRASG12D was significantly attenuated by mutation at K104. Biochemical and structural analysis indicated that the G12D and K104Q mutations cooperate to suppress GEF-mediated nucleotide exchange, explaining the preferential effect of K104Q on oncogenic KRAS. Furthermore, K104 functioned in an allosteric network with M72, R73, and G75 on the α2 helix of the switch-II region. Intriguingly, point mutation of glycine 75 to alanine (G75A) also showed a strong negative regulatory effect on KRASG12D. These data demonstrate that lysine at position 104 is critical for the full oncogenic activity of mutant KRAS and suggest that modulating the sites in its allosteric network may provide a unique therapeutic approach in cancers expressing mutant KRAS. SIGNIFICANCE: An allosteric network formed by interaction between lysine 104 and residues in the switch-II domain is required for KRAS oncogenicity, which could be exploited for developing inhibitors of the activated oncoprotein.

RAS-dependent RAF-MAPK hyperactivation by pathogenic RIT1 is a therapeutic target in Noonan syndrome-associated cardiac hypertrophy.

RIT1 is a RAS guanosine triphosphatase (GTPase) that regulates different aspects of signal transduction and is mutated in lung cancer, leukemia, and in the germline of individuals with Noonan syndrome. Pathogenic RIT1 proteins promote mitogen-activated protein kinase (MAPK) hyperactivation; however, this mechanism remains poorly understood. Here, we show that RAF kinases are direct effectors of membrane-bound mutant RIT1 necessary for MAPK activation. We identify critical residues in RIT1 that facilitate interaction with membrane lipids and show that these are necessary for association with RAF kinases and MAPK activation. Although mutant RIT1 binds to RAF kinases directly, it fails to activate MAPK signaling in the absence of classical RAS proteins. Consistent with aberrant RAF/MAPK activation as a driver of disease, we show that pathway inhibition alleviates cardiac hypertrophy in a mouse model of RIT1 mutant Noonan syndrome. These data shed light on the function of pathogenic RIT1 and identify avenues for therapeutic intervention.

Machine Learning-Driven Multiscale Modeling: Bridging the Scales with a Next-Generation Simulation Infrastructure.

Interdependence across time and length scales is common in biology, where atomic interactions can impact larger-scale phenomenon. Such dependence is especially true for a well-known cancer signaling pathway, where the membrane-bound RAS protein binds an effector protein called RAF. To capture the driving forces that bring RAS and RAF (represented as two domains, RBD and CRD) together on the plasma membrane, simulations with the ability to calculate atomic detail while having long time and large length- scales are needed. The Multiscale Machine-Learned Modeling Infrastructure (MuMMI) is able to resolve RAS/RAF protein-membrane interactions that identify specific lipid-protein fingerprints that enhance protein orientations viable for effector binding. MuMMI is a fully automated, ensemble-based multiscale approach connecting three resolution scales: (1) the coarsest scale is a continuum model able to simulate milliseconds of time for a 1 µm2 membrane, (2) the middle scale is a coarse-grained (CG) Martini bead model to explore protein-lipid interactions, and (3) the finest scale is an all-atom (AA) model capturing specific interactions between lipids and proteins. MuMMI dynamically couples adjacent scales in a pairwise manner using machine learning (ML). The dynamic coupling allows for better sampling of the refined scale from the adjacent coarse scale (forward) and on-the-fly feedback to improve the fidelity of the coarser scale from the adjacent refined scale (backward). MuMMI operates efficiently at any scale, from a few compute nodes to the largest supercomputers in the world, and is generalizable to simulate different systems. As computing resources continue to increase and multiscale methods continue to advance, fully automated multiscale simulations (like MuMMI) will be commonly used to address complex science questions.

GAP positions catalytic H-Ras residue Q61 for GTP hydrolysis in molecular dynamics simulations, complicating chemical rescue of Ras deactivation.

Functional interaction of Ras signaling proteins with upstream, negative regulatory GTPase activating proteins (GAPs) represents a crucial step in cellular decision making related to growth and survival. Key components of the catalytic transition state for Ras deactivation by GAP-accelerated hydrolysis of Ras-bound guanosine triphosphate (GTP) are thought to include an arginine residue from the GAP (the arginine finger), a glutamine residue from Ras (Q61), and a water molecule that is likely coordinated by Q61 to engage in nucleophilic attack on GTP. Here, we use in-vitro fluorescence experiments to show that 0.1-100 mM concentrations of free arginine, imidazole, and other small nitrogenous molecule fail to accelerate GTP hydrolysis, even in the presence of the catalytic domain of a mutant GAP lacking its arginine finger (R1276A NF1). This result is surprising given that imidazole can chemically rescue enzyme activity in arginine-to-alanine mutant protein tyrosine kinases (PTKs) that share many active site components with Ras/GAP complexes. Complementary all-atom molecular dynamics (MD) simulations reveal that an arginine finger GAP mutant still functions to enhance Ras Q61-GTP interaction, though less extensively than wild-type GAP. This increased Q61-GTP proximity may promote more frequent fluctuations into configurations that enable GTP hydrolysis as a component of the mechanism by which GAPs accelerate Ras deactivation in the face of arginine finger mutations. The failure of small molecule analogs of arginine to chemically rescue catalytic deactivation of Ras is consistent with the idea that the influence of the GAP goes beyond the simple provision of its arginine finger. However, the failure of chemical rescue in the presence of R1276A NF1 suggests that the GAPs arginine finger is either unsusceptible to rescue due to exquisite positioning or that it is involved in complex multivalent interactions. Therefore, in the context of oncogenic Ras proteins with mutations at codons 12 or 13 that inhibit arginine finger penetration toward GTP, drug-based chemical rescue of GTP hydrolysis may have bifunctional chemical/geometric requirements that are more difficult to satisfy than those that result from arginine-to-alanine mutations in other enzymes for which chemical rescue has been demonstrated.

RASopathy mutations provide functional insight into the BRAF cysteine-rich domain and reveal the importance of autoinhibition in BRAF regulation.

BRAF is frequently mutated in human cancer and the RASopathy syndromes, with RASopathy mutations often observed in the cysteine-rich domain (CRD). Although the CRD participates in phosphatidylserine (PS) binding, the RAS-RAF interaction, and RAF autoinhibition, the impact of these activities on RAF function in normal and disease states is not well characterized. Here, we analyze a panel of CRD mutations and show that they increase BRAF activity by relieving autoinhibition and/or enhancing PS binding, with relief of autoinhibition being the major factor determining mutation severity. Further, we show that CRD-mediated autoinhibition prevents the constitutive plasma membrane localization of BRAF that causes increased RAS-dependent and RAS-independent function. Comparison of the BRAF- and CRAF-CRDs also indicates that the BRAF-CRD is a stronger mediator of autoinhibition and PS binding, and given the increased catalytic activity of BRAF, our studies reveal a more critical role for CRD-mediated autoinhibition in BRAF regulation.

Structure of the SHOC2-MRAS-PP1C complex provides insights into RAF activation and Noonan syndrome.

SHOC2 acts as a strong synthetic lethal interactor with MEK inhibitors in multiple KRAS cancer cell lines. SHOC2 forms a heterotrimeric complex with MRAS and PP1C that is essential for regulating RAF and MAPK-pathway activation by dephosphorylating a specific phosphoserine on RAF kinases. Here we present the high-resolution crystal structure of the SHOC2-MRAS-PP1C (SMP) complex and apo-SHOC2. Our structures reveal that SHOC2, MRAS, and PP1C form a stable ternary complex in which all three proteins synergistically interact with each other. Our results show that dephosphorylation of RAF substrates by PP1C is enhanced upon interacting with SHOC2 and MRAS. The SMP complex forms only when MRAS is in an active state and is dependent on SHOC2 functioning as a scaffolding protein in the complex by bringing PP1C and MRAS together. Our results provide structural insights into the role of the SMP complex in RAF activation and how mutations found in Noonan syndrome enhance complex formation, and reveal new avenues for therapeutic interventions.

Exploring CRD mobility during RAS/RAF engagement at the membrane.

During the activation of mitogen-activated protein kinase (MAPK) signaling, the RAS-binding domain (RBD) and cysteine-rich domain (CRD) of RAF bind to active RAS at the plasma membrane. The orientation of RAS at the membrane may be critical for formation of the RAS-RBDCRD complex and subsequent signaling. To explore how RAS membrane orientation relates to the protein dynamics within the RAS-RBDCRD complex, we perform multiscale coarse-grained and all-atom molecular dynamics (MD) simulations of KRAS4b bound to the RBD and CRD domains of RAF-1, both in solution and anchored to a model plasma membrane. Solution MD simulations describe dynamic KRAS4b-CRD conformations, suggesting that the CRD has sufficient flexibility in this environment to substantially change its binding interface with KRAS4b. In contrast, when the ternary complex is anchored to the membrane, the mobility of the CRD relative to KRAS4b is restricted, resulting in fewer distinct KRAS4b-CRD conformations. These simulations implicate membrane orientations of the ternary complex that are consistent with NMR measurements. While a crystal structure-like conformation is observed in both solution and membrane simulations, a particular intermolecular rearrangement of the ternary complex is observed only when it is anchored to the membrane. This configuration emerges when the CRD hydrophobic loops are inserted into the membrane and helices α3-5 of KRAS4b are solvent exposed. This membrane-specific configuration is stabilized by KRAS4b-CRD contacts that are not observed in the crystal structure. These results suggest modulatory interplay between the CRD and plasma membrane that correlate with RAS/RAF complex structure and dynamics, and potentially influence subsequent steps in the activation of MAPK signaling.

Asynchronous Reciprocal Coupling of Martini 2.2 Coarse-Grained and CHARMM36 All-Atom Simulations in an Automated Multiscale Framework.

The appeal of multiscale modeling approaches is predicated on the promise of combinatorial synergy. However, this promise can only be realized when distinct scales are combined with reciprocal consistency. Here, we consider multiscale molecular dynamics (MD) simulations that combine the accuracy and macromolecular flexibility accessible to fixed-charge all-atom (AA) representations with the sampling speed accessible to reductive, coarse-grained (CG) representations. AA-to-CG conversions are relatively straightforward because deterministic routines with unique outcomes are achievable. Conversely, CG-to-AA conversions have many solutions due to a surge in the number of degrees of freedom. While automated tools for biomolecular CG-to-AA transformation exist, we find that one popular option, called Backward, is prone to stochastic failure and the AA models that it does generate frequently have compromised protein structure and incorrect stereochemistry. Although these shortcomings can likely be circumvented by human intervention in isolated instances, automated multiscale coupling requires reliable and robust scale conversion. Here, we detail an extension to Multiscale Machine-learned Modeling Infrastructure (MuMMI), including an improved CG-to-AA conversion tool called sinceCG. This tool is reliable (∼98% weakly correlated repeat success rate), automatable (no unrecoverable hangs), and yields AA models that generally preserve protein secondary structure and maintain correct stereochemistry. We describe how the MuMMI framework identifies CG system configurations of interest, converts them to AA representations, and simulates them at the AA scale while on-the-fly analyses provide feedback to update CG parameters. Application to systems containing the peripheral membrane protein RAS and proximal components of RAF kinase on complex eight-component lipid bilayers with ∼1.5 million atoms is discussed in the context of MuMMI.